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Abstract 

New exactly solvable problems have already been studied by using a modification 
of the factorization method introduced by Mielnik. We review this method and 
its connection with the traditional factorization method. The survey includes the 
discussion on a generalization of the factorization energies used in the traditional 
Infeld and Hull method. 



1 Introduction 

Since the begining of the Quantum Mechanics (QM) there exists only a narrow set 
of exactly solvable physical problems. This set includes the well known harmonic 
oscillator and hydrogen-like potentials, the Morse potential, the square well potential 
and a few others. The study of such problems has been traditionally done by using 
the method of orthogonal polynomials and the factorization method. In principle, it 
is no matter the method used in the construction of the analytical solutions to the 
Schrodinger equation, however, the factorization (introduced by Schrodinger ^ and 
Dirac [Q]) avoids the use of cumbersome mathematical tools and it has been succes- 
fully applied as an elegant approach in solving essentially any problem for wich there 
exists exact solution. This is clear from the classification and study of the factorizable 
second order differential equations given by Infeld and Hull [g]. It is remarkable that, 
although the factorization method was developed to solve the eigenvalue problem re- 
lated with the time-independent Schrodinger equation, it is also a very powerful tool 
in the study of the recurrence relations obeyed by special functions m. 
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Until the begining of the 1980's decade, it was a very common opinion that the 
method was completely explored. The change of mentality became after the semi- 
nal work of Mielnik [Q], The fundamental idea introduced by Mielnik in such paper 
is to consider not the particular but the general solution to the Riccati type equa- 
tion connected with the Schrodinger-Dirac-Infeld-Hull (SDIH) approach. Hence, the 
Mielnik's method takes into account an important term of the factorization opera- 
tors unnoticed by the traditional SDIH method. The first application of this modified 
method was developed in the study of the harmonic oscillator by Mielnik himself, 
giving rise to the construction of a one-parameter family of new exactly solvable po- 
tentials, different than the harmonic oscillator potential but having the same spectrum 
as it. Early on, the Mielnik method was applied by Fernandez in the construction 
of a one-parameter family of new exactly solvable radial potentials, isospectral to 
the hydrogen-like radial one [|[|. In the same year Nieto gave the links between this 
modified SDIH method and the supersymmetric (SUSY) approach to QM [Q]. Its con- 
nection with Darboux transformation [|]] has been discussed by Bagrov and Samsonov 
[pO (see also QKJ), and some links with the inverse scattering method have been given 



by Zakhariev et al. [ 1 1 1. By now the Mielnik method has been maturing for one and a 
half decade, and far from becoming a closed chapter in solving the Schrodinger equa- 
tion, it is continuously reviewed and applied in diverse prospects [12-16]. For more 
historical links and related approaches see the interesting account given by Rosu in 
this volume. 

In this contribution we shall present a general survey of the factorization method. 
Our main interest is to state clearly the differences between the SDIH and the Mielnik 
methods and, at the same time, to give the introductory material allowing their ap- 
plication in the derivation and study of new exactly solvable problems. In Section 2 
we will review the SDIH method. Section 3 deals with the Mielnik method. Finally, 
in Section 4, we shall review the further generalizations of the Mielnik's approach 
recently introduced in [ |l4] , |l5| ]. Each section contains the corresponding harmonic 
oscillator and hydrogen-like potentials as particular cases. 



2 Traditional factorization 



The Schrodinger equation for a time-independent one-dimensional system is given by 
the usual eigenvalue equation 

s2 



ft 
2m 



dx 2 



v(x) 



ip{x) = £ip(x) 



(1) 



where the potential v(x) can be a singular (e.g. the hydrogen-like potential) as well as 
a nonsingular (like the free particle or the harmonic oscillator potentials) real func- 
tion, and the wave function ip{x) must be square-integrable, that is, the integral 
/ | ip(x) | 2 dx is finite. One often uses wave functions which are normalized, such 
that: 

r \i>(x)\ 2 dx = l, (2) 
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equation (g) is usually called the normalization condition. 

As is well known, any physical problem must state not only the differential equa- 
tion which is to be solved but also the boundary conditions which the solution must 
satisfy. Moreover, it is remarkable to notice that the determination of solutions sat- 
isfying of the boundary conditions is often as difficult a task as the solving of the 
differential equation. At the present state, it is apparent the essential role played by 
the normalization condition (g) in quantum theory because it represents the strongest 
boundary condition to be satysfied by the solutions to (jj]). 

The spirit of the factorization method, as introduced by Schrodinger and Dirac, 
is to write the second order differential operator H in ([l]) as the product of two first 
order differential operators a and a', plus a constant e. The form of such operators 
depends explicitly on the particular potential v(x) we are dealing with, and implicitly 
on the factorization energy e, as we shall see below. We are going to sketch now a 
brief survey of the SDIH method. In order to eliminate a cumbersome notation let us 
rewrite the Schrodinger equation (fy) in the dimensionless form 



dx 2 



V{x) 



^(x) = E1>(x), (3) 



where V(x) = 2mv(x)/h 2 and E = 2m£/h 2 . The factorization method aims to find 
two operators 

a=—+f3 p (x), a} = --^+/3 p (x), (4) 
ax ax 



such that the Hamiltonian 



can be written as 



2 



H=~ + V(x) (5) 
ax z 



H = a t a + e. (6) 

It is a matter of substitution to show that (3 p (x) has to satisfy the following Riccati 
equation: 

-*3' p (x)+(3l(x) = V(x)-e. (7) 

The SDIH method takes into account only the particular solution to the Riccati 
equation (0), this is the reason why we have labeled the /3-functions, from the very 
begining, with the subindex "p". In the following table we reproduce the Infeld-Hull 
results for the harmonic oscillator as well as for the hydrogen-like potentials. The 
traditional representation for the harmonic oscillator is recovered from (Q) and (||) by 
take H H/2, a -> (1/V2)a, and a f (1/V2)a f . 
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Potential 


Factorization 
energy 


Particular solution to 
the Riccati's equation 


V{x) = x 2 


e=l 

"modified method" [|[ 


(3 p (x,e) = x 




for a fixed Z 


P P (r,e)= l ~-\ 
for a fixed I 



Let us remark that, in general, given any e 6 R for which there exists a solution 
to ([7]), there is a factorized expression for H as it is given in However, it is 
worthwhile noticing that, because the Riccati equation (^) is a non-linear differential 
equation, it is not, in general, integrable by quadratures. It therefore defines a family 
of transcendental functions [|l7|, and the search of an adequate factorization energy 
permitting the construction of the /3 p -functions is not an easy task. In fact, Infeld and 
Hull have reported only one factorization energy for each potential V(x), hence, they 
have given the first step in the systematic construction of exactly solvable problems in 
Quantum Mechanics. 

Now, let us take the product (||) in the opposite order and, by using (Q) and ([7]), we 

get 

oat + e = H + 2/3' p (x). (8) 

In a first view, one might be tempted to belive that equation (||) brings out a new 
Hamiltonian H + 2(3' p (x), but it is not always the case, as we are going to see in the 
next sections. 



2.1 The harmonic oscillator revisited 

As is well known, the solutions to (||), with V(x) — x 2 , are given by Hermite poly- 
nomials type wave functions ip n (x), whose corresponding eigenvalues are given by 
E n = 2n + 1, n = 0, 1, 2, ... The Infeld-Hull results are, in this case, e = 1 and 
(i p {x) = x, hence we have a — d/dx + x, and equation ^ can be rewritten as 

d 2 

H=-—,+x 2 = aa)-l. (9) 
dx 

It must be clear that, in this case, equation (^|) does not give any new Hamiltonian. 
Now, by using (||) and the corresponding factorization we get the well known 
commutation rules satisfied by the operators a, and H: 



[a, a 1 ] = 2, [H, a] = -2a, [H, a f ] = 2a f . 



(10) 



On the factorization method in Quantum Mechanics 



289 



Furthermore, it is well known that, for the harmonic oscillator problem, the operators 
of factorization and a are also the corresponding raising and lowering operators for 
the eigenfunctions of the Hamiltonian. 

On the other hand, the method by itself gives us the chance of constructing a first 
solution to (|^). In order to make clear this argument let us consider the equation (^), 
and suppose that there exists a function iPn(x), such that aipisr(x) = 0, hence we 
will have Hi/jn(x) — eipN(x). The solution to this first order differential equation, 
providing the solution to (0) as given, is 

^(i)«e-r*W (11) 

If iPn{x) satisfies the normalization condition (||), then it is a wave function with 
eigenvalue e, else it is not a physical solution to the Schrodinger equation ([}]). For the 
system we are dealing with in this section, the function iPn(x) = Co exp(— x 2 /2) is 
a physical eigenfunction of H with the eigenvalue e = En = 1, hence N = 0, and 
4>n(x) = ipo{x) is the ground state wave function of the harmonic oscillator. 

Let us consider now the Hamiltonian (^). The solution to the first order differential 
equation a'tpM (%) — 0, given by iPm(x) oc exp(x 2 /2), is an eigenfunction of H 
with eigenvalue Em = E-\ = — 1 ^ {E n }. It is clear that iPm(x) is an unphysical 
solution because either it not satisfies the normalization condition (||), or because its 
eigenvalue Em is not allowed in the spectrum of H. We want to remark again that 
these unphysical solutions are usually discarded by the traditional methods of solution 
in Quantum Mechanics. In Section 4 we will have opportunity to show their relevance 
in the construction of new exactly solvable potentials. 



2.2 The hydrogen-like potential revisited 

The standard procedure to deal with hydrogen-like potentials in QM reduces to solve 
the eigenproblem for a particle in a one-dimensional effective potential Vj(r) = 1(1 + 
l)/r 2 — 2/r, where I = 0, 1, 2, is the azimuthal quantum number and r is a dimen- 
sionless radial coordinate. By simplicity, instead of working with the standard radial 
wavefunctions R l n (r), we will work with the functions ip n .i(r) = rR l n (r), with an 

inner product defined by (ip n ,h' l Pn,l') = ^ n Jo + °° , /w( r )V'n,J'( r ') < ^ r ' < 00 ■ As it is 
well known, the eigenvalues of the radial Hamiltonian 

ar z r z r ar z 

for a fixed I, are given by 



E n =E lM =- — , N2 ; fc = l,2,3,... (13) 



1 

(I + k) 

where I + k = n. Now, for the sake of clarity, let us take the following notation for 
the Infeld-Hull factorization operators of the Hamiltonian (|l2j) 

d I 1 + d I 1 

ai = — + - - y, a] = -— + - - -, (14) 

dr r I dr r I 
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hence, the corresponding factorizations (|6|) and (|8|) can be written as 

a\ai-j2=Hu aiaj - = iZj-i. (15) 

Notice that the right hand side equation (|5j) (e.g. (||)), can be rewritten as 

, 1 

a w a i+i ~ (; + 1)2 = Hi. (16) 

In the present case, the operators of factorization a \ and a;, map the solutions of the 
corresponding equation ([|), with a given energy, into solutions to (|j) with the same 
energy but changing the value of the azimuthal quantum number I. Hence, a\ and a/ 
are respectively the raising and lowering operators for the quantum number I, besides 
we have i> n ,i+i{r) oc aj +1 ip n> i(r), and ip n ,i-i(r) oc ajip n j(r). 

Let us conclude this section with the observation that the physically relevant first 
solution to (|j), derived from the SDIH method for the Hamiltonian ([j~2]), is given 
by ij}M (Q oc r l+1 exp(— r/(l + 1)), which is solution of a| +1 -0j\,/(r) = 0. Hence, 
from @ and (]l3|), it is apparent that it is an eigenfunction of Hi with eigenvalue 
e = Eis = —1/(1 + l) 2 . On the other hand, the unphysical solution is now given by 
V>jv(r) oc r~ l exp(r/(7)), with the forbidden eigenvalue E^o = —l/l 2 . 



3 Mielnik's factorization 

In this section we shall review the modification of the SDIH method as introduced by 
Mielnik. As we have discussed in the introduction, the Mielnik's method aims to find 
two operators 

A=-^+P(x), At =- A +0(x), (17) 
ax ax 

factorizing the Hamiltonian (|j) in the form 

H = A i A + e, (18) 
where the function j3(x) is the general solution of the Riccati equation 

-P'(x) +(3 2 (x) = V(x) - e. (19) 

It is well known that, given any particular solution (3 v (x) to (^9[), the corresponding 
general solution is given by two successive quadratures [n7|, and one can write 



0(x) =(3 p {x)--^ky\\- I e 2 f^ d ydx 



(20) 



where A is an integration constant. 

Now, by taking the factorization ( |l8| ) in the opposite order we get 

Atf +e = H + 2(3' p (x)-2-^\n^\- J e 2 f W . (21) 
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From Sections 2.1 and 2.2, it is clear that only the term H + 2/3' (x), in the right 
hand side of (j2l|), would correspond to the initial Hamiltonian H. The element of 
arbitrarness becomes then from the logarithmic term in such equation. Hence, we can 
define a new Hamiltonian H = ~(d 2 /dx 2 ) + V(x), where 

V[x) = V{x) + 20' (x) (22) 

and such that 

H = AA ] + e. (23) 



Equations ( 22 ■ 23|) open a new interesting problem: Is the following eigenvalue equa- 
tion exactly solvable? 

H$(x) = E4>{x). (24) 

The answer is affirmative and the corresponding solutions can be constructed by using 
the solutions of (||), with H given in (|l8|), and the following theorem: 

Intertwining Theorem 

Let ip( x ) be an eigenfunction of H with eigenvalue E, then Aip(x) ^ is an 
eigenfunction of H with eigenvalue E. Similarly, if <j)[x) is an eigenfunction of 
H with eigenvalue E, then A^<j)(x) ^ is an eigenfunction of H with eigen- 
value E. 
Proof. 

H[Aip(x)] = A[A^A + e]ip(x) = E[A^(x)\, 
H[A^(x)\ = A ] [AA ] + e]<j>(x) = E[A ] 4>{x)]. Q.E.D. 

Then we have ip(x) oc Aip(x), and ip(x) oc A^4>(x). Notice that, in the present 
approach, the operators of factorization A? and A, are not the corresponding raising 
and lowering operators for the eigenfunctions of either H or H. Now, as we have 
discussed in Section 2, if the solutions to (Q) are wave functions, then they satisfy the 
normalization condition (||). Hence, we have 

(if, ft) = (Atfj, Atp') = (A^Aii, ft) = ((H - e)^, ft) = {E - e)(V>, ft) 

and, if ijj(x) £ L 2 (R), we can take ip(x) — (E — f)~ x / 2 Aip(x). It is interesting to 
notice that, although we have used all the eigenfunctions of < ]18| ) in the construction 
of the new set {ip(x) = (E — e)~ 1 / 2 Aip(x) \ ip(x) £ L 2 (R)}, it is not, in general, a 
complete set. Suppose there exists a function tpMe (%) sucn that it is orthogonal to all 
tp(x), then we will have 

$ Me ,$) = (E-ey^i^AiP) = (E- e)-^ 2 (A^ Me ,^) = 0. 

Therefore, the solution to the first order differential equation A^ipMeix) = 0, rep- 
resents a possible missing element in the new basis {ip(x)}. Such solution is given 
by 

feWoce/^W*, (25) 



292 



J.O. Rosas-Ortiz 



and by using (J20j) we can rewrite it as 

VMe{x) oc TWTTTj = 1 12 ( \J ' (26) 

A-/VJ ^ dz dy X ~I 1>m(v)*V 

where ^m{x) is solution of a'ipM(x) = 0. 

Now, if ( p5| ) satisfies the normalization condition (^), then it is a physical solution 
to the eigenvalue problem (p4j), with eigenvalue e (see equation (p3|)), and it has to 
be added to the new basis {ip(x)}. In general, the eigenfunctions of H are given 
by {ipMe(x)} U {i/j(x)}. Moreover, the Hamiltonian H is almost isospectral to H, 
because its spectrum is the same as the spectrum of H plus a new level at e: {E} = 
{e, -E 1 }. In other words, we have solved the eigenvalue equation ( |24| ) by using the 



intertwining theorem and by considering the missing state ( |25[ ). 

Let us remark that the eigenproblem connected with H can be first solved by using 
the SDIH method and after by the Mielnik approach, which gives no new solutions for 
H, but allows the construction of new exactly solvable Hamiltonians H. Therefore, 
the Mielnik method is a further step in the systematic construction of exactly solvable 
problems in QM. 



3.1 First order intertwining approach 

As we have dicussed in the last section, the Mielnik's approach allows one to derive 
new Hamiltonians H departing from the factorization of a given Hamiltonian H. It 
has been also shown that the eigenfunctions of H are an adequate point of departure in 
solving the eigenproblem connected with H. Now, a natural question arises: Are the 
Hamiltonians H and H connected by a more direct relation? The answer is afirmative 
and it is clearly stated by the algebraic tool called the first order intertwining method 
[|l8p. Now, we are going to show how the basic results of the intertwining method can 
be recovered from the Mielnik's approach. With this aim, let us apply the operator A 
to the right of equation (^3|), hence, we get 

HA= [AA* +e]A = A[A i A + e] = AH. (27) 

Equation ( p7[ ) means that the Hamiltonians H and H are connected by the operator 
A. Notice that the multiplying order of the operators plays a fundamental role. Sup- 
pose now that H is a Hamiltonian with known solutions to the corresponding equation 
, then ( p7[ ) immediately leads to the first part of the intertwining theorem. In the in- 
tertwining/argon it is said that A is an intertwiner operator because it intertwines (the 
eigenfunctions of) the operators H and H. In this case the Hamiltonian H is called 
the initial one, while H is the intertwined Hamiltonian, which is to be determined by 
the method itself. On the other hand, the operator A^ intertwines the Hamiltonians H 
and H in the form: 



HA^ = [A^A + e]A^ = A^[AA^ +e] = A^H. 



(28) 
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In this case, considering the multiplying order of the operators, it is suitable to take H 
as the initial Hamiltonian, therefore, equation ([^8]) allows one to recover the second 
part of the intertwining theorem. On the other hand, it is a matter of substitution to 
show that equations ( |l9| ) and (22) are sufficient conditions validating the intertwining 



relationship (E7h. Moreover, all these equations lead, in a natural form, to the factor- 



ization of the initial Hamiltonian H with the multiplying order of (18). In a similar 
form, the sufficient conditions to satisfy ( ^8| ) are given by 

f3'{x)+ (3 2 {x) = V(x)-e (29) 

and 

V(x) = V(x)-2j3'(x), (30) 

permitting the factorization of the Hamiltonian H as it has been made in (|23"|). 

In the last section we have done the assumptions validating the intertwining rela- 
tionship (]27[). Now, we shall proceed in the opposite way, i.e., we will suppose that 
the eigenvalue problem ( p4| ) has been solved by the SDIH method and, therefore, we 
will focus on solving the corresponding equation for H . With this aim, let us consider 
now the function a(x) as a particular solution to (p^), then the corresponding general 
solution is given by 

(3{x) = a(x) + -^ln^y + J e~ 2 f a{v)dv dx^ , (31) 

where 7 is an integration constant (compare (|T]) with (j20|)). The corresponding eigen- 
functions are given by those derived from the intertwining theorem, plus the missing 
state 

i> Ne (x) oce-J ™*> = . (32) 

. rx -2 1 a(z)dz , 

7 + J e J v ' ay 

Let us remark that, although it is no matter what kind of factorization (e.g. equa- 
tions (18]) or ([0])) we use for the initial Hamiltonian in the intertwining approach 



(|27[|28), it is usually important to choose it adequately, in order to get a physically 
permissible intertwined Hamiltonian. Last statements will be clear in the next subsec- 
tion, where we are going to apply the Mielnik method to the harmonic oscillator and 
to the hydrogen-like potentials by using the intertwining approach. 



3.1.1 Classical isospectral potentials 

Harmonic oscillator 

We take the initial Hamiltonian H = —d 2 /dx 2 + x 2 + 2 as it has been factorized in 
(8|^). We want to factorize it in the form (p3|), hence, by making V(x) = x 2 + 2 in 
(29|), we get the particular solution a(x) = x (see the table of Section 2). Therefore, 
the corresponding intertwined potential ( j30| ) is given by 
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If I7I > y/n/2, the above potential has no singularity and behaves like x 2 for x — > 
±00; and so, one obtains here a one-parameter family of self-adjoint Hamiltonians H 
in L 2 (R) The first eigenfunctions of H are given by the intertwining theorem, 
while the corresponding missing state ( |32| ) is given by 

e -x 2 /2 

ll>Ne(x) OC -j—, (34) 

7 + Jo e v d V 

which is a square-integrable function because its behavior for x — ► ±00, provided that 

| 7 | > 

Hydrogen-like potential 

Let us take iJ; as it has been factorized in ([j^). We are going to factorize it in the form 
(|l8|). By substituting Vi(r) in ( |l9| ) we get the particular solution /3 p (r) = //r — 1 /I, 
with fixed I (see the table in Section 2), and the corresponding intertwined potential 
Vi-i(r) (see equation (E2)) is given by 




If A; > (2/)! (l/2) 2L+1 , or A ; < 0, for a fixed I, the third term has no singularities. 
Furthermore, for r — > +00, we have VJ_i(r) = V;_i(r), and therefore we obtain 
a one-parameter family of self-adjoint Hamiltonians Hi-\ in i 2 (R) [^|. The first 
eigenfunctions of H1-1 are given by the intertwining theorem, while the correspond- 



ing missing state (26 ) is given by 



« Xt-xluU/'dy' (36) 

which is a square-integrable function because its behavior at r — 0, and in the limit 
r -> +00, provided that A; > (2/)! (Z/2) 2 ' +1 , or A; < 0. 



4 Further factorizations 

In the last sections we have done a survey of the factorization methods available for 
the construction of analytical solutions to the Schrodinger equation. Now, we remark 
that, although the Mielnik method is a powerful tool in the derivation of new Hamilto- 
nians whose corresponding eigenproblem is analytically solvable, not all the particular 
solutions to the Riccati equation reported in |^] have been taken into account for their 
generalization through the Mielnik's approach. On the other hand, it is a matter of 
fact that the Infeld-Hull results in the table of Section 2 (e.g. |j^]) correspond to only 
one factorization energy e for each potential. Therefore, there is only one SDIH fac- 
torization available for the corresponding Hamiltonian and we have at hand only one 
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family of isospectral potentials connected with it. In the present section we are going 
to sketch some steps allowing the selection of new factorization energies. With this 
aim, we shall take H in (JTsb as the initial Hamiltonian [n9|. The transformation 



(3{x) 



dx 



hiu(x), 



leads ( |l9| ) to the second order differential equation 



dx 2 



V{x) 



u(x) = eu(x). 



(37) 



(38) 



Hence, the /3-functions are directly connected with the eigenfunctions of the initial 
Hamiltonian H, corresponding to the eigenvalue e. Notice that, if the factorization en- 
ergy e belongs to the spectrum {E}, then the solutions u(x) to ( |38| ) are the physically 
relevant functions ip(x). Now, from Section 3.1.1, we know that the factorization en- 
ergy e = 1, for the harmonic oscillator H + 2, does not belong to E n+ i. On the other 
hand, for the hydrogen-like potential we have e = — l/l 2 , which is not allowed in 
Ei k, K > 0. Hence, following [^[ and [^], we shall consider just the cases where the 
factorization energy e does not belong to the spectrum of H. Therefore, the solutions 
to ( |38| ) do not have direct physical meaning, but we have just seen that they natu- 
rally lead to the factorization of the Hamiltonians H and H in the spirit of Mielnik's 
approach, providing as well the explicit form for the new potentials V(x). 

In order to solve ([38|) we take the factorization of H as given in I 
two options: 

(0 a u(x) = 0, 

(if) au(x) = w(x) 0; a^au{x) = a) w(x) = 0. 



, hence we have 



Notice that u(x) = iPn(x), with iI)n{ x ) given in (|1 1[), is the solution of (/), whereas 
w(x) = iPm(x) — Ci/iPn(x), with c\ a constant, is the solution to the right hand side 
equation (if). Hence, the general solution to the left hand side equation (if) is given by 



u(x) = 4> N {x) 



Co + ci 



(39) 



where cq is an integration constant. Now, by using (37) and <\l 1|), with c\ = — 1 and 
Co = A, we get (|o|). The general solution of ( |l9| ) is then constructed by the general 
solution of (Bq), with the unphysical eigenvalue e. 



4.1 New isospectral potentials 



In this section we are going to solve, explicitly, the eigenvalue problem (38) for the 
harmonic oscillator and the hydrogen-like potentials. Hence, our results will bring out 
new factorization energies permitting the construction of new one-parameter families 
of potentials isospectral to each one of the potentials mentioned before. 
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4.1.1 The hydrogen-like potential 



In order to solve (|38|), with V(x) = Vj(r), let us write the factorization energy as 
e[ fe ' = —1/(1 + k) 2 , k ^ K, I > 0. It is clear that, for a fixed value of I, we have 
e| fe ' 7^ V fc 7^ K. Now we make the transformation 

u | fc )(r) =r-'e r /( z+fc )$ ; (fe) (r), (40) 

leading to a confluent hypergeometric equation for \r), whose general solution 
for the discrete values k = 0,-1, —2, — (I — 1) is given by the linear combination 
of confluent hypergeometric functions [E0| 



$[ fc) (r) = 1 F 1 [fe,-2/,-2r/(/ + fc)] 

- n fc [2r/(i + fc)] 1+2 'iFi[l + k + 21, 2 + 21, -2r/(i + k)], 



with 



_ r(i + ifcj) r(-2Q w 
r(2? + 2) r(-2i + |*p ' ' 



The general solution to the corresponding Riccati equation ( |i9| ) arises after introduc- 
ing © in (|7]), which gives: 



^w-j-jn^-s^M. <>»■ 

The intertwined potentials V^^r) have the same singularity at r = as V;_i(r), 
provided that A; takes only values in the domain [0|: 



(—oo,l), for \k\ even; 
(1, 00), for |fc| odd. 



By taking k = and aJ "* = 0, the last results reduce to the corresponding SDIH re- 
sults for the same value of k and Aj '' = (2/)!(Z/2) 2/+1 7 / ~ 1 , we get the Fernandez 
results [^] and, finally, the corresponding Abraham and Moses results [|2|] arise for 
k = and A^ 0) -> 1, with 2 = 1. 



4.1.2 The harmonic oscillator 

For the harmonic oscillator potential we take V(x) = x 2 as the initial potential, hence, 
in order to solve the corresponding equation (|3S|), let us make the transformation 



u(x) = cxp(-x 2 /2), 



(41) 
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and by introducing instead of x the variable y = x 2 we then find $(?/) to satisfy a 
confluent hypergeometric type differential equation, whose general solution is given 



by [20| 



m = by) + " ^tS f 1/! W T' bv)- 

v 4 J 



The general solution to the corresponding Riccati equation ( |29| ) arises after the trans- 
formation f3(x) = (d/dx) \nu{x), which gives 

/3(x) = -x- -f-ln$(z 2 ). 

dx 

The intertwined potential becomes free of singularities if e < 1 and \u\ < 1 [ |l4[ pl| ]. 
In particular, we can take = — 2fc - 1, with k = 0,1, 2, ...[Q. 

In this case, the SDIH results arise by taking k = and i/ = 0, the Mielnik results 
[Q| are recovered with k = and z^ * 1 = (y / 7r/2)7 _1 and, by taking the same value of 
k, with v — ¥ 1, the Abraham and Moses family is recovered. We conclude this Section 
remarking that the factorization energies e\ k \ and e^ k \ generalize the choice made for 
the SDIH and Mielnik factorizations of the hydrogen-like and harmonic oscillators 
potentials respectively. 



4.2 Higher order factorizations 

The Mielnik factorization allows us to construct further new families of exactly solv- 
able potentials departing from any exactly solvable problem and iterating after the 
procedure. The new families will inherit the number of parameters from the old ones. 
The keystone of that construction is given by [114, 0, Hffl: 



I3 n+1 {x) = -0 n (x) - [ — ^- ) , n = 1,2, . 

,Pn+l — Pn 



where n denotes the number of iterations. 

Notice that the last equation works immediately for e n +x ^ e n . Hence, by giv- 
ing n different factorization energies one can iterate n times the Mielnik's approach 
constructing nth-parametric families of analytically solvable potentials. The results 



presented in this paper reproduce those recently derived in [14 |15|. Although we 
have shown how this method works by solving the discrete problems connected whith 
two interesting physical systems, its applications include the study of systems with 
continuous spectrum (e.g. e„+i — > e„) 
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